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Abstract 

This paper presents the geometric setting of quantum variational principles and extends it to 
comprise the interaction between classical and quantum degrees of freedom. Euler-Poincare reduc¬ 
tion theory is applied to the Schrddinger, Heisenberg and Wigner-Moyal dynamics of pure states. 
This construction leads to new variational principles for the description of mixed quantum states. 
The corresponding momentum map properties are presented as they arise from the underlying uni¬ 
tary symmetries. Finally, certain semidirect-product group structures are shown to produce new 
variational principles for Dirac’s interaction picture and the equations of hybrid classical-quantum 
dynamics. 
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1 Introduction 


After Kibble’s investigation [27] of the geometric properties of quantum state spaces, geometric for¬ 
mulations of quantum dynamics have been attracting much attention over the last decades [3, 4, 6, 7, 
8, 13, 14, 15, 18, 21, 34, 39, 26, 42]. In turn, geometric quantum dynamics has opened several modern 
perspectives: for example, Fubini-Study geodesics have been introduced in Grover’s quantum search 
algorithms [33, 45] and in time-optimal quantum control [10, 11, 16], while the holonomy features 
arising from the quantum geometric phase [1, 7] and its non-Abelian extensions [2, 14] have been 
proposed in quantum computation algorithms [22, 31, 40]. 

Most approaches deal with pure quantum states and involve the geometry of the Hopf fibration 

'Ip I—>■ , 

where S (denotes the unit sphere in a complex Hilbert space Jif (so that ip € S (^) is a unit vector 
in J^), while is the corresponding projective space containing the projections := ip'ipK The 
geometry of the above Hopf bundle is well known and has been widely studied in the finite dimensional 
case = C”, although some studies extend to consider inhnite-dimensional Hilbert spaces [13]. In 
the hnite dimensional case, one can emphasize the symmetry properties of the Hopf bundle by writing 

S{C^) =U{n)/U{n-l), PC^ = S{C^)/U{1) = U{n)/{U{n - 1) x U{1)) . (1) 

The interplay between the geometry of the Hopf bundle and its symmetries is the basis of geometric 
quantum dynamics. For example, the emergence of principal bundles leads to the usual horizontal- 
vertical decomposition in terms of a principal connection that is strictly related to Berry’s geometric 
phase. This is a beautiful picture, whose symplectic Hamiltonian properties have been widely inves¬ 
tigated after Kibble’s work [27]. 

In this paper, we aim to present how this geometric framework emerges naturally from the unitary 
symmetry properties of quantum variational principles. Time-dependent variational approaches have 
been most successful in chemical physics (here, we recall the celebrated Car-Parrinello model in 
molecular dynamics [9]). The most fundamental quantum variational principle is probably due to 
Dirac and Frenkel (DF) [17, 19]. This action principle produces Schrodinger equation ihip = Hip as 
the Euler-Lagrange equation associated to 

6 f (ip,itnp — Hip^ dt = , 

Jti 

where H is the quantum Hamiltonian operator and we introduce the pairing {A, B) and the inner 
product 

{A\B) = Tr{A^B) , so that (A, B) := Re {A\B) 

and Im {A\B) = {iA,B). Various properties of the variational principle above have been studied over 
the decades [30, 28, 36], after it was first proposed in the context of Hartree-Fock mean field theories. 
For example, it is known that the DF action principle is simply the quantum correspondent of the 
classical Hamilton’s principle on phase space <5 [p ■ q — H{q, p)) dt = 0, so that h{ip, idip) acquires 
the meaning of canonical one form on (analogously, p • dq is the canonical one form in classical 
mechanics). However, to our knowledge, an investigation of the geometric symmetry properties of 
quantum variational principles has not been carried out. Although an early attempt was proposed in 
[28], the emergence of the Hopf bundle in this context has not been presented so far. For example, 
the momentum maps associated to quantum variational principles have never been considered in 
the literature, while they are essential geometric features often associated to fundamental physical 
quantities. Even in the simplest situation, the phase invariance of quantum Lagrangians produces the 
momentum map identifying the total quantum probability HV'lp. 
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The present work applies well known techniques in the theory of geometric mechanics [32], which 
focus on the symmetry properties of the dynamics. For example, momentum map structures are 
seen to emerge: some are new, while others are related to the principal connections associated to 
quantum geometric phases. Within geometric mechanics, we shall be using the specihc tool of Euler- 
Poincare theory [24] that typically applies to variational principles with symmetry. As we shall 
see, besides recovering well known relations, these theory allows us to formulate new variational 
principles for various quantum descriptions, such as the Liouville-Von Neumann equation, Heisenberg 
dynamics, Moyal-Wigner formulation on phase space and the Ehrenfest theorem for the evolution of 
expectation values. Some of these descriptions of quantum mechanics have been lacking a variational 
structure, which is now provided in this paper for the first time. Here, we shall not dwell upon various 
complications that may emerge in infinite dimensional Hilbert spaces M’ and we assume convergence 
where necessary. When convenient, we shall consider dynamics on finite dimensional spaces and rely 
on the possibility of extending the results to the infinite dimensional case. 


2 Euler-Poincare variational principles in the Schrodinger picture 

This section presents the Euler-Poincare formulation of quantum dynamics in the Schrodinger picture. 
Two main examples are considered: the Schrodinger equation as it arises from the Dirac-Frenkel theory 
and the Fubini-Study geodesics. Their geometric features will be analyzed in terms of momentum 
maps. 


2.1 Euler-Poincare reduction for pure quantum states 

Upon denoting by the tangent bundle of the Hilbert space consider a generic Lagrangian 

L = L{'ilj,tp) , (2) 

so that the assumption of quantum evolution restricts ip to evolve under the action of that unitary 
group V({J^), that is 

m = U{t)Pjo, U{t) € (3) 

where pJo is some initial condition, whose normalization is ordinarily chosen such that JIV’D IP = 1- 
Then, V’o S implies ^p{t) € at all times. 

The relation (3) takes the Lagrangian to a Lagrangian of the type Lf), which then 

produces Euler-Lagrange equations for the Lagrangian coordinate U G Moreover, by following 

Euler-Poincare theory [24], one denotes by u(J^) the Lie algebra of skew Hermitian operators and 
defines 

i{t) ■.= U{t)U-^{t) Gu(Jf’). 

Since ip = one obtains the reduced Lagrangian 

£ : u{J£f) xJ£f^R, P;) := L(V’, p'lp) 


and the Euler-Poincare variational principle 


^ 'ip)dt = 0 . 

Jti 


Then, upon computing 


6p = r] + [rj, V] , S^p = rjip 


where r] := {5U)U one obtains the following result. 


( 4 ) 

( 5 ) 
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Theorem 1 Consider the variational principle (4) with the auxiliary equation ip = and the varia¬ 
tions (5), where rj is arbitrary and vanishes at the endpoints. This variational principle is equivalent 
to the equations of motion 


_d ^ _ r ^ 

TtSf ~ 


1 

2 


6ip 


ip"^ — Ip 




( 6 ) 


Here, we use the ordinary definition of variational derivative 

SF(<,) 

for any function(al) F £ C°°{M) on the manifold M. In typical situations, the reduced Lagrangian 
is quadratic in ip, so that the ZY(l)-invariance under phase transformations takes the dynamics to the 
projective space Indeed, as we shall see, the reduced Lagrangian ip) can be written typically 

in terms of the projection = ipip"^ to produce a new Lagrangian 


I : u{J^) X ^ R, 

In this case, a direct calculation shows that 

^Pip \'niPij )\) pp [CjP'i/'] 
and the previous theorem specializes as follows 


Theorem 2 Consider the variational principle 5p^) dt = 0 with the relations (7) 
r) + [r],(,], where rj is arbitrary and vanishes at the endpoints. This variational principle is 
to the equations of motion 


and 6^ = 
equivalent 


d 61 

■ SF 


■ 61 

dt Sf, 

6^ 

— 

Spp^'^ _ 


Pp K) Pp] ■ 


( 8 ) 


Then, the unitary symmetry properties of the Lagrangian naturally take the evolution to the correct 
quantum state space (for pure states), that is the projective space FM’. In the following sections, 
we shall specialize this construction to two particular examples and we shall present the momentum 
map properties of the underlying geometry as well as their relation to the usual principal connections 
appearing in the literature. 


2.2 Dirac-Frenkel variational principle 

It is easy to see that upon following the construction from the previous section, the DF Lagrangian 

L{ip,ip) = (^ip,ihip — Hip) (9) 

produces the Euler-Poincare variational principle 

6 f {ip,ihf^ip — Hip) dt = 0 

Jti 

For simplicity, here we are considering a time-independent Hamiltonian operator H. Then, upon 
computing 


^ = 2{ih£, - H)ip, ^ = -ihipip "^, 
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the first of (6) yields 


[(ihe-//),V'V’^] = o. (10) 

Upon setting = C” and making use of the anticommutator bracket {A,B} = AB + BA, the 
solution of the above equation can be written as 

^ + ih~^H = {1 — ( 11 ) 

for an arbitrary time-dependent skew Hermitian matrix K{t). The meaning of this solution will be 
clear in Section 2.5. In the Schrodinger picture, the relation (11) recovers the usual phase arbitrariness, 
as it is shown by simply using the second in (6) to write 

ihtl) = H^|J + a'lf ), ( 12 ) 

where a{t) := 2h{ip^,K) (notice that we have chosen a unit initial vector so that H'fAolP = = !)• 

In the above equation, the term a'lp generates an arbitrary phase factor. Then, equation (12) can be 
easily written in the form of a projective Schrodinger equation [27] 

(1 — — Hip) = 0 , 

which can be recovered from the constrained DF Lagrangian [36] 

L{ip,ip,X,X) = {ip,ihip - Hip) + AdlV'll^ - 1). (13) 

As we shall see, the right hand side of (11) modifies the usual Heisenberg picture dynamics. 

2.3 Mixed states dynamics and its Wigner-Moyal formulation 

It is easy to see that all the phase terms in the previous section are consistently projected out by simply 
defining the Lagrangian l{p,,pp) = {pp,ih^ — H) so that the hrst of (8) reads [{ihS, — H),p^] = 0 and 
the second recovers the quantum Liouville equation for pure states. 

In the remainder of this section we shall generalise the previous approach to consider a new 
variational principle for mixed quantum states. Let us consider the Lagrangian 

l{i,p) = {p,iK-H) (14) 

where ^ = UU~^ as before, while p is a density matrix undergoing unitary evolution p(t) = UpoU^■ 
In the case of mixed states, we have p^ ^ p although the trace invariants Tr(p”) are still preserved. 
We notice that a simple computation yields 

=Vih]: P= [?, P], (15) 

and therefore the application of Euler-Poincare theory is straightforward. Then, one obtains precisely 
the same equations as in (8) (upon replacing pp by p), which in turn give 

H),p]=0, p=[^,p]. (16) 

At this point the Liouville-Von Neumann equation 

ihp = [H, p] 

is obtained by direct substitution. Notice that the solution of the first equation in (16) differs from 
(11), since p^ ^ p. For example, one has particular solutions of the form ih^ — H = 

This reflects the very different geometric structures underlying mixed states and pure states. For a 
geometric description of mixed states in terms of coadjoint orbits and orthogonal frame bundles, we 
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refer the reader to [34]. We emphasize that the action principle associated to the Lagrangian (14) is 
very different from the one proposed in [23] and to our knowledge it has not appeared before. 

Motivated by applications in chemical physics, we show how the above variational principle recovers 
the celebrated Wigner-Moyal picture of quantum dynamics on phase space. This formulation [35, 43] 
is based on the Weyl correspondence between linear operators and phase space functions (see e.g. 
[44]). For simplicity, this section presents the Euler-Poincare formulation on the two-dimensional 
phase space (one spatial dimension), however this can be easily generalized to higher dimensions. 

Consider an arbitrary linear operator A € the corresponding phase-space function is given 

by the Wigner transform a{x,p) = W(M) and the latter can be inverted by using the Weyl transform, 
A = W“^(a). More explicitly, one has 

yV{A):= — dx'{x + x'\A\x — x')e 
TTh J 

= 2 / dxdx' \x + x'){x — x' 


dp a{x,p)e « 


Then, the Moyal bracket (see [35, 44] for its explicit definition) is defined in such a way that the 
commutator between two quantum operators is taken into the Moyal bracket of the corresponding 
phase space functions, that is [44] 

{{a,6}) = l>V([A-Bl). 

At this point, one can express the Lagrangian (14) in terms of phase space functions. Indeed, upon 
defining the Wigner distribution W{x,p) = Vy(/o) and by replacing the inverse relation p = 'yV~^{W) 
in (14) one obtains the equivalent variational principle on phase space 


<5 



dxdp W{x,p) [hT{x,p) — H{x,p)J dt = 0 , 


(17) 


where we have defined H{x,p) = W{H) and T{x,p) := W(z^). Then, upon recalling the relations 
(15) and the first in (5), one computes 

(iT = ^ + n{{e, T}}, sw = h{{e, w}}, dtW = h{{T, w}}, 

where 0 := W{iri), and the Euler-Poincare variational principle (17) gives 


dW{x,p,t) 

di 


{{H{x,p),W{x,p,t)}} . 


Again, we notice that the variational principle (17) has never appeared before in the literature and 
it is very different from other approaches proposed earlier, such as [37]. In particular, the variational 
principle (17) is entirely derived from the Dirac-Prenkel Lagrangian and no assumption has been made 
other than unitary evolution. 


2.4 Geodesics on the space of quantum states 

Eubini-Study geodesics on PC"' are used in various situations of quantum mechanics. Their applica¬ 
tions in quantum search algorithms [33, 45], time-optimal control problems [10, 11, 16] and holonomic 
quantum computation [22, 31, 40] emphasizes their importance and makes this example especially 
interesting. In the general case, geodesics are optimal curves in the sense that they minimize the 
distance between two quantum states. 
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The Fubini-Study geodesics are defined as geodesic equations on minimizing the the Fubini- 
Study distance. These geodesic flows can be written explicitly as Euler-Lagrange equations associated 
to the action principle d dt = 0 with Lagrangian L : TJif —>■ M given by (see e.g. [18]) 


2 \m^ 


(18) 


More explicitly, lengthy computations yield 

(llV'f - (IIV’ll^'0 - 2(V’l'i/’)V’) = 0. 


Notice that this approach does not involve normalized vectors in Jff. However, a simple way to recover 
normalization is to use a constrained Lagrangian such as that in (13). 

Application of the Euler-Poincare theory to Fubini-Study geodesics can be performed again by 
following the procedure outlined in Section 2, without modifications. Then, upon recalling that is 
Hermitian, one obtains the reduced Lagrangian (set h= 1, for convenience) 

KC, P^) = -^ ((pi/., ■ (19) 

Notice that we recover the well known relation l{^,Ptp) = 1/2 ((T^) — (T)^) [1], where we have intro¬ 
duced the energy operator £ = and we have used the standard expectation value notation. Then, 
upon replacing the variational derivatives 

^ = “2 ^ = 2^PP^^} - {pplOPi’ ( 20 ) 

in equations (8), standard matrix computations give 

^ Pip'j 0 5 P^ll [^) P^}] 5 

where the first emphasizes the following conservation form of the Fubini-Study geodesic equation 

^ ((1 - 2p^)p^ = ^ = 0 , 

as it arises from the left-invariance of the Lagrangian (18) (see e.g. [5]). Notice that applying the 
above conservation law to and writing ('0|V') = i{i'ip,'ip) yields 

v) - [tjj -b 

and since (V’lV’) = “('l/IV’)! expanding leads to 

(l — {pp — 2{'il)\'ip)ip) = 0 . 


This geodesic flow can be also recovered as an Euler-Lagrange equation by adding a normalisation 
constraint AdlV'jp — 1) to the Lagrangian (18) (cf. the constrained DF Lagrangian (13)). 

Geodesic flows on the quantum state space have always raised questions concerning their underlying 
geometric properties [3, 42]. For example, in holonomic quantum computing, a fundamental role is 
played by the connection form {'ip\'ip), whose loop integral defines the celebrated geometric phase. In 
addition, connection forms also allow the usual horizontal/vertical decomposition on the Hopf bundle. 
In geometric mechanics, this decomposition can be performed by using a more sophisticated theory 
than Euler-Poincare reduction. This is called Lagrange-Poincare reduction [12] and it was recently 
formulated in the context of homogeneous spaces (arising from symmetry breaking) in [20] . Without 
entering the technicalities of Lagrange-Poincare reduction, we shall only mention that this theory often 
takes advantage of a particular connection form that appears to have a precise physical meaning in 
many different cases: this is called the mechanical connection and it is defined in terms of a momentum 
map, another fundamental object in geometric mechanics. This is the topic of the next section. 


7 



2.5 Momentum maps of quantum variational principles 

The first momentum map one encounters in quantum mechanics is probably the density matrix for 
pure states [15]. More precisely, the action of the unitary group hl{J^) on the quantum Hilbert space 
^ (endowed with the symplectic form = 2^ (zV'i, " 02 )) produces the momentum map 

J('0) = G , 

as it can be easy obtained by the general formula (J('0),^) = [32, 25], holding for an 

arbitrary skew-hermitian operator ^ G u(J^). Also, restricting to consider phase transformations 
yields the total probability or, more precisely, the quantity J('0) = 

Other than those above, other momentum map structures appear in geometric quantum dynamics 
and each correspond to different group actions and different reduction processes. It turns out that in 
quantum variational principles, the most important momentum map is associated to the action of the 
isotropy subgroup of the initial state. In order to explain this statement, let us replace the relation 
(3) in a Lagrangian of the type (2) and observe that this produces a Lagrangian L : TU{J^) ^ M by 

J:{U,U) := LiU^PoMo) ■ (21) 

Although, this Lagrangian is not symmetric under right multiplication, i.e. 

C{U,U)^ C{UU',UU'), U' 

the invariance property is recovered by restricting to the isotropy group of tpo, that is 

= {U^ UiM’) I UiPo = V’o} • 


Indeed, one evidently has 


C{U,U) = £{UUo,UUo), ( 22 ) 

and one may choose the initial vector ipQ to coincide with the basis vector = (0 • • • 0 1)^^. 

Remark 1 (Analogy with the heavy top dynamics) We observe that the argument above holds 
in a wide range of situations, including, for example, the Lagrangian reduction for the heavy top 
dynamics [20]. In that context, the unitary group is replaced by the rotation group 50(3) and the 
isotropy symmetry is defined to preserve the gravity veetor, thereby leading to planar rotations in 
50(2). The Noether’s conserved quantity (i.e. the momentum map) is then the vertical angular 
momentum. 

At this point, it is natural to ask what the momentum map is for the reduction of quantum 
variational principles. More particularly, we look for the momentum map associated to (the cotangent 
lift of) the right action of on the cotangent bundle T*U{J^). In the general case, it has 

recently been shown [20] that the momentum map for the right action of a subgroup Go C G on (the 
trivialisation of) the cotangent bundle T*G ~ G x g* reads 

J(g,;x) = G(Ad» (23) 

where {g,p) G G x g*. Ad*// = g^iJ.g~'^ is the standard matrix coadjoint representation and l* is the 
dual of the Lie algebra inclusion / : go g. To simplify the treatment, set = C"" and choose 
-00 = (0 ... 01)^ without loss of generality. Then, U,Pq{J^) = U{n — 1) C U{n) = hl{M’) and the 
group inclusion Z//(n — 1) hl{n) is 


U[n — 1) 3 Uq 


Uo 0 
0 1 


G U{n). 


(24) 


The corresponding Lie algebra inclusion t : u(n — 1) u(n) reads 


^(■Co) = 


0 


0 

0 


while its dual i* is given by r*(//) = (1 — /9^g)/i(l — that is the standard projection on the 

upper left block. This result is independent of the number of dimensions and it leads to the following 
momentum map: 




1 

2 


{(1 — 2p^g), Ad^ p] + I Ad^ p)ppq 


= Ad^ 





(25) 


where we have simply rewritten (23) by replacing the formula for i*. Here, we recall ip^^ = Ad^(ipp), 
from the definition pp := = Uppf^U~^. Therefore, because of the symmetry property (22) 

possessed by any Lagrangian of the type (2), the corresponding Euler-Poincare equations (6) conserve 
Ji{U,5i/6^). More particularly, one shows that any quantum system with an arbitrary Lagrangian of 
the type (21) produces dynamics on the zero-level set of Ji. This is easily shown by using the relation 
= L('0,^'i/’) =• ■^(^iV’)) so that 




y Sip ^4’ ) 


(26) 


thereby verifying Ji{U,6£/6^) = 0. 

Now that we have characterised the momentum map associated to the action of ^/(n — 1), we recall 
that all physically relevant Lagrangians must be also phase invariant, so that they can be eventually 
written in terms of the projection p^ = G Therefore, the most general symmetry group 

of the Lagrangian (21) has to include phase transformations and this leads us to consider the direct 
product U{n — 1) x hl{l). The latter can be embedded in U{n) by the inclusion 


U{n — 1) xU{l) 3 {Uq, (/?) !->• 


Uo 0 
0 


£U{n), 


(27) 


where the minus sign in the exponent is of purely conventional nature. Since the momentum map 
associated to U{n — 1) has already been presented in (25), we need to compute only the momentum 
map associated to the group 7/(1) C 7/(n), endowed with the group inclusion 

Z^(l) [I ^ 


Upon computing the dual of the corresponding Lie algebra inclusion i(a) = —iap^^ € u(n) and by 
identifying u(l) ~ M, one has the momentum map formula 

MU, p) = i {ppo I Ad^ p) . (29) 

Any quantum system with an arbitrary Lagrangian of the type (21) takes the above momentum map 
to the form 

as it is easily shown by using the relation (26). Here, the momentum map 77 : T*S'(C”) —)■ M arises 
from the action of 7/(1) (notice that we idenitified u(l) ~ M) on the cotangent bundle T*S{C^). For 
example, the DF Lagrangian yields {ip, 6L/d'lp) = while one verifies that the FS Lagrangian 

(18) makes 77 vanish identically. 
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In more generality, the momentum map corresponding to the action of the full symmetry group 
U{n — 1) X U{\) is given by 

= {Ji{U,n),J2{U,n)) = Q{(1 - 2/9v>o),Ad^/z} + {p^JAdu i {pp^ \ Ad^/r)^ (30) 

Therefore, because of the symmetry property (22) possessed by any phase-invariant Lagrangian of the 
type (2), the corresponding Euler-Poincare quantum dynamics (8) conserves the momentum map J 2 
and lies on the kernel of Ji. 

Remark 2 (Phases and W(l)—actions) Notice that this section has used a U{1) —action that is 
different from usual phase transformations. Indeed, while the latter act on vectors by the diagonal 
action if i->- the U{1) —action used in this section reads if ^ {1 — -|- e“®‘^ppg)V'- However, 

usual phase transformations (denoted by Wd(l) to emphasise the diagonal action) are a subgroup of 
b({n — 1) X U{1), as it is given by the inclusion 

Uffl) 9 (^) ^ ( ) G W(n - 1) X U{1) C U{n ). 

Therefore, our treatment naturally includes the ordinary phase transformations, whose correspond¬ 
ing momentum map is given by i*ffJi{U,p),J 2 {U,p)). Here, is the dual of the inclusion : 
Urf(l) ^ u(n — 1) X u(l) given by id{oi) = (—zq:1,q:) (again, we identify Urf(l) ~ M). Upon using 
the pairing {{p,uj),{p,a)) = {p,p) + ua, a direct calculation shows L*^{p,uj) = uj-\-Tr{ip), so that 
L*{jffU,p),J2{U,p))=iTT{p). 

The momentum maps provided in this Section are of paramount importance in geometric quan¬ 
tum dynamics, as they incorporate essential geometric properties. For example, it is interesting 
to notice that the momentum maps (25) and (30) can be used to define connection forms, respec¬ 
tively on the bundles U{n) —>■ S'(C"') and U{n) —>■ PC”, as they are given in (1). These con¬ 
nection forms are obtained by applying Lagrange-Poincare reduction for symmetry breaking and 
we refer the reader to [20] for more details on this topic. In the particular case under consid¬ 
eration, one identifies u(A:)* ~ u(/c) so that the dual of the inclusion U : u(n)* —)• u(n — 1)* de¬ 
termines a projection P : u(n) —> u(n — 1), which in turn can be used to define the mechani¬ 
cal connection A{U) = P(Ad^-i ^) = (1 — ppg)(Ad[/-i ^)(1 — ppg) on the bundle U{n) —>• S'(C”). 
An analogous construction yields a connection on the bundle lA{n) —)■ PC”, given by A{U) = 
(1 — ppg)(Adf/-i ^)(1 — ppg) -|- i(Ad[ 7 -i f I ppo)ppo. Also, it is well known that A{if) = {if,iif) is a 
principal connection on the Hopf bundle 5(C”) ^ PC”. Then, the momentum maps presented in 
this section generate a connection form on each bundle of the diagram below. The study of these 
connection forms and their curvatures is left for future work. 


U{n) 



The arguments in this Section clarify the meaning of equation (11). Indeed, the latter can be 
interpreted as simply saying that the infinitesimal generator ^ G u(n) (or, equivalently, the Hamiltonian 
operator H) is dehned only up to an element of the isotropy subalgebra u(n — 1) x u(I). This can be 
made explicit upon introducing (kq, oq) G u(n— 1) xu(l) such that {1 —2/9.^, «:} = (! — p^)ko(1 “ Pip)^- 
iaoPi/'- Although computing if = ^if returns the usual phase arbitrariness, as shown in equation (12), 
the relation (11) discloses a rich geometry content underlying quantum evolution. Indeed, it reminds 
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that the propagator of quantum dynamics is defined up to elements of Z^(n — 1) x U{1) (not only up 
to phases in which is actually a non Abelian symmetry group (unlike phase transformations). 

As we shall see in the next Section, this argument accounts for propagators that depend on the initial 
quantum state ipQ. 

3 The Heisenberg and Dirac pictures 

While the previous sections mainly dealt with the Schrodinger picture of quantum mechanics, the 
Heisenberg picture is rather unexplored in the geometry of quantum evolution. For example, one is 
interested in the role of projection operators, as they emerge from the projective geometry of the 
quantum state space in the Schrodinger picture. We shall consider the Heisenberg picture for the DF 
Lagrangian and the Fubini-Study geodesics. 

3.1 Euler-Poincare reduction in the Heisenberg picture 

Most of this Section is devoted to the Heisenberg picture for the DF Lagrangian. It is easy to see that 
(9) can be written in the Heisenberg picture by introducing 

in ■= U-^U = ^, Hh:=U^HU. 

Indeed, with these definitions, the Lagrangian (9) becomes 

K^h,Hh) = {P' 4 , 0 , - Hh) ■ (31) 

Then, upon computing 

= PH - [pH,iH], 5 Hh = [Hh,Ph] , Hh = [Hh,^h] (32) 

(with r]H := U~^5U), inserting the Lagrangian (31) in the variational principle Hh) dt = 0 

yields the following Euler-Poincare equations: 

[ih^H - Hh, P' 4 , 0 ] = 0 • (33) 

At this point, we observe that although the above relation is satisfied by ih^n — Hh = otl (so that 
the Heisenberg Hamiltonian Hh is defined up to a phase factor a), more general solutions are present 
such as 

= -ifH^Hn + {1 - 2/3^0, k} (34) 

(k being arbitrary and skew Hermitian, see equation (11)). These solutions have the property of 
depending on the initial state tpo. Upon setting Jif = for simplicity, one may choose = 
(0 .. .01)^ without loss of generality. Interestingly enough, these more general solutions lead to the 
unfamiliar equation 

Hh = [HH,{d-‘iP4,o,K'}] , (35) 

so that the Hamiltonian operator Hh is not conserved in the general case. Although this may seem 
surprising, we observe that the above dynamics does not change the physics of the system under 
consideration. For example, we observe that the total energy is preserved: 

{Hh) = {p4,o\Hh) = 0, 

as shown by a direct verification. Moreover, one realizes that the above dynamics of quantum Hamil¬ 
tonians returns exactly the Schrodinger equation (12): indeed, one has 

= -ih~^HH'ipo - 2(/0^ol^)V’o 
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so that, recalling if; = Uipo and applying U on both sides returns (12). Notice that it is indeed essential 
that identifies the initial quantum state. We conclude that the physical content is unaltered by 
the Heisenberg equation (35), which in turn generalizes the standard Heisenberg dynamics (recovered 
by At = 0) to incorporate the geometry of quantum dynamics. 

As a practical example, we consider spin dynamics in the Heisenberg picture. In this case, the 
Hamiltonian reads Hh = n • Sh, where Snit) = U{t)~^SU{t) in standard spin operator notation. 
The DF Lagrangian (9) is written in the Heisenberg picture as 1{Ch, Sh) = — n • Sh), so 

that the Euler-Poincare equations 

[ih^H - n • Sh, Ppo] = 0, ith = [-S'h, C/a] 


specialize to yield 

Sh = ih~^[n ■ Sh + /{2ppo - 1, At},^^] =nx Sh - [{1- 2ppg, At}, Sh] 


Notice that this approach can be applied in the general case. For example, one can study linear 
oscillator dynamics by recalling the Hamiltonian Hh = huja^j^an and following precisely the same 
steps as above. The present approach leads to the following Heisenberg equation: 

Ah = ih~^ [Hh, Ah] - [{1 - , At}, Ah] , 

where At is an arbitrary skew-symmetric operator and Hh undergoes its own evolution (35). In 
addition, equation (34) yields a new form of the propagator equation 

U = ih~^HU + U{1 — 2ppq,kI} . (36) 


It is necessary to point out that, since At is an arbitrary skew-Hermitian matrix parameter, one can 
simply choose it in such a way that {1 — 2/?^g,At} = 0, thereby eliminating the dependence of the 
propagator on the initial conditions. A similar argument leads to eliminating phase terms in the 
Schrodinger equation (12) [27]. 

In all this Section, we assumed the initial state is a pure state ijjQ. If this is not the case, then 
different solutions of the type C// = —ih~^HH + ac (with [ac, po] = 0) are allowed by equation (33), 
because in this case pp^ is replaced by a density matrix po 7^ Po- 

The Heisenberg picture is particularly natural for the description of FS geodesics. Indeed, while 
in the Schrodinger picture the right unitary symmetry is broken by pp, in the Heisenberg picture one 
may use the full left symmetry of the Lagrangian (18), that is Upon 

recalling (3) and by setting h = 1 for convenience, one obtains the Euler-Poincare Lagrangian 

A^h) = --^(^{pPoI^h) + {ppoAinf^ ■ (37) 


Then, upon using the first of (7), the variational principle K^h) dt = 0 yields 


{C/AjPPo} “ ‘^{ppo\iH)Ppo + 2(ppolC//)C//,Ppo 


= 0 , 


which reflects all the properties already discussed in Section 2.4 


3.2 Dirac-Frenkel Lagrangian in the Dirac picture 

This section extends the arguments from the previous Sections to formulate a new variational prin¬ 
ciple for quantum dynamics in the Dirac (interaction) picture. As we shall see, the Euler-Poincare 
construction involves the semidirect product of the unitary group with itself. 
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In the Dirac picture, the Schrodinger Hamiltonian operator is split in two parts as H = Hq + Hi, 
where Hq is typically a simple linear Hamiltonian, while Hi usually contains nonlinear potential terms. 
This replacement can then be inserted in the DF Lagrangian (9). However, it is convenient to keep 
track of the quantum state tps that is propagated by Hq, such that ihips = Hqi/js (up to phase terms). 
Then, one is led to consider the following DF Lagrangian 

- {Hq + Hi)'ijj) + {'ijjs,ihiijs - Hoips)- (38) 

Performing Euler-Poincare reduction by replacing the evolution relation (3) on the first part yields 
the Lagrangian 

- Hq- Hi) + - HqiPs) , 

with = U and ^ := UU~^. At this point, the propagator associated to Hq can be used 
to replace the evolution relation tpsit) = in the second term, thereby leading to the Euler- 

Poincare Lagrangian 

Pp^Hqj, Hij) = {ppj,ih^i - Hqj - Hij) -h {pp^,ih^Q - Hqj) (39) 

where p^^ = V’oV'o nnd we have introduced the following definitions 

V-/ = Do-V, Hjj = U^^H,Uq, (q = Uo^Uq, e/ = C/o”'et/o = Ad^-ie. (40) 


In order to write the resulting equations of motion, we start by using the last two definitions in 
(40) to compute the variations 

^(^ 0 , 6 ) = (% + [Co,'no],'ni + [^0,'ni] - [vo,(i] + [ni,^i]) (4i) 

where rjQ = Uq^6Uq and rp = Ad^-l{{5U)U~^). One recognizes that the variations (41) are Euler- 

Poincare variations of the type Su = C + [Cj where [•, -Jg is the Lie bracket on g = uq{J^)(^u{J^), 
that is the Lie algebra of the semidirect product group lJQ{Jlf)®U{Jlf). Here, the group UQ^Jlf) 
is a copy of the unitary group and one thinks of Uq{J^) as accounting for the propagators 

Do GLfo(^). 

Computation of the other variations by using (40) yields 


5ppi = [pi - PV'/] ’ m ]) 


so that the variational principle 1 {^q,^j, pj, Hqj, Hij) dt = 0 produces the following equations of 
motion for an arbitrary Lagrangian t. 


d 61 
'dtJij ~ 

\f “1 

+ 


+ 

61 ■ 

Ppll r 

L oppj j 

= 0, 

d 61 

1 

+ 


- 

>0 Q. 

•Cl 

C!. 

1 _ 1 

- 

\h 

[ 6Hqp 


Pi = [^/ - Co, Pi] , Hjj = [Hjj, ^o] . 



61 

J^i 


= 0 , 


(42) 

(43) 

(44) 


Then, computing the variational derivatives of the Lagrangian (39) and replacing them into (42) and 
(43) gives 

Ppi] = [Ha,I + Hij, ppj] , [ih^Q - Hqj, p^^] = 0 
As seen in Section 3.1, the second relation above is solved by 

^0 = -ifr^HQj -h {1 - k} , 

so the second in (44) gives 

Hq,i = [Hq,i, {1 - 2p^^,k}] , Hip = ih~\HQp, Hip] -|- [Hip, {1 - 2p^^,K}] . 
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On the other hand, the first in (44) becomes 

= ihr^ [pp,, Hij] + [pp^, {1 - 2p^p, k}] . 

Then, we notice that the choice k = 0 returns the usual quantum dynamics in the Dirac picture. As 
we know from Section 3.1, the K-terms do not change the overall physical content of the dynamics. 
For example, a direct calculation verifies the following energy conservations: 

^ + Hij) = 0 , ^ (p^g \Ho,i)=0. 

As we have seen, the application of Euler-Poincare reduction theory reveals the geometric features 
emerging in the Heisenberg picture of quantum dynamics. These geometric features reveal the form 
(36) of the propagator equation, without affecting the physical content of quantum dynamics. When 
this form of the propagator equation is considered in the Dirac picture, this introduces extra terms in 
the dynamics, which still preserve the total energy of the system. 

In the next Section, we shall formulate a new variational principle for the interaction of classical 
and quantum degrees of freedom. 


4 Classical-quantum variational principles 

The interplay of quantum and classical degrees of freedom has always attracted much attention in 
quantum mechanics. For example, the consistent formulation of hybrid quantum-classical models 
in molecular dynamics remains an outstanding issue [38]. In this section, we present a geomet¬ 
ric formulation of the most elementary system coupling classical and quantum dynamics. This is 
given by combining the Ehrenfest equations for the expectation of the canonical variables with the 
Schrodinger/Liouville equation for the quantum degrees of freedom. More particularly, we shall present 
a novel variational principle for the Ehrenfest mean field model and in more generality for expectation 
value dynamics. 

4.1 The classical-quantum mean field model 

In order to approach the dynamics of quantum expectations, we observe that the mean field closure 
of any classical-quantum system can be derived in first instance by the following Lagrangian 

L(z, z, '0, '0) = ^ z • Jz -b (V’, ihip - H (z)'0), (45) 

where —Sijdz^Adz^ is the canonical symplectic form and H{z) is a Hermitian operator depending on 
the classical degrees of freedom z = (x, p). Indeed, the corresponding Euler-Lagrange equations yield 

z = J (V’|Vzll(z)'0) , ih'il; = , (46) 

thereby recovering the ordinary mean field model of classical-quantum dynamics (see e.g. equations 
(12.2)-(12.4) in [41]). Here, purely classical dynamics is recovered by the phase type Hamiltonian 
H{z) = h(z)/ (here, I denotes the identity operator on while purely quantum dynamics is 

recovered when Vzi4(z) = 0. 

An Euler-Poincare formulation of the above equations can again be obtained by letting the quan¬ 
tum state evolve under unitary transformations. This leads to a coupled Euler-Lagrange equation for 
(z, z) and the Euler-Poincare equations for the quantum dynamics (expressed in terms of either "0 or 
its density matrix). However, in order to find a full set of Euler-Poincare equations that includes the 
classical evolution, we may choose to evolve the phase-space vector z under the action of the Heisen¬ 
berg group (i.e., phase-space translations), which is prominent in the theory of quantum coherent 
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states. To this purpose, consider a curve h{t) = (h(t), ^{t)) in the Heisenberg group ~ 

and let the phase space vector z evolve as 


z(t) = Zo + h(t). 


(47) 


Also, we recall ip(t) = Then, upon inserting the auxiliary phase factor (p in the Lagrangian 

(45), the latter becomes 

Lzo,-,Po{h,h,U,U) = ^ h • JI(zo + h) + {U'lpo, {ihU + pU - H{zo + h)t/)V'o) • 

The above Lagrangian is of the type 

L^o,^o ■ TniR^^) X M 


and its dynamics can be approached by Euler-Poincare reduction. Therefore, in order to find an 
expression for the reduced variable C := hh~^, we dehne the Lie algebra element 


iCA) ■= 


ds 


{g{s)h € 


2n\ 


s=0 


where g{s) = (g(s),i9(^)) S ?7(M^"') is a curve such that ^(0) = h and 5 ^( 0 ) = h (for some fixed time). 
Here, we recall the Heisenberg group operation 


gh = l^g + h, ?? + + - g • Jh 

which gives h~^ = (—h, —<p). Eventually, one finds 


Vg,h€n(R^A, 


(48) 


C={CA)= ( h, (/i - • Jh ) , 


so that the Euler-Poincare Lagrangian is written as 


^iCAA, 2 ,p^) = C • Jz + {pp,if>^ + 4>- H{z)) , (49) 

where we have used the convenient initial condition zq = 0 in (47), without loss of generality. Notice, 
this Lagrangian is of the type 

i : (f)(M2^) © u(Jf')) X X PJf) ^ M . 

Then, the Euler-Poincare equations follow in the theorem below, upon recalling the infinitesimal 
adjoint representation 

'(' 2 ) = (0,—Cl ■ JC 2 ) (50) 

in the Heisenberg Lie algebra f)(]R^”'). The present treatment is now extended to the case of mixed 
quantum states. 


Theorem 3 Consider the variational principle 


6 ^C • Jz + (p, zhC + 0 — L7(z))^ dt = 0 


and the variations 

((5C,5C>) = ( 7 , 0 + C-J 7 ) > = ^z = 7 , 6p = [n,p], 

where ( 7 , 0 ) and g are arbitrary and vanish at the endpoints. Together with the auxiliary equations 

z = C , P=[^,P], 

this variational principle is equivalent to the equations of motion 

z = J(p|VziL(z)), ihp=[H{z),p] . 
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Proof. Consider the general Lagrangian of the form p^). By direct substitution of the 

variations into the variational principle 


ft2 

Jti 


K 


, 7 ) + 




M 

Sz 


, 7 ) + 


Tp' 


p]))dt = 0, 


one writes the Euler-Poincare equations as 

dt 6C S(j) Sz ’ dt d(j) 
In particular, for the Lagrangian (49), we have 
Si 6£ Si 

= JZ, — = (p|l), 


= 0 . 


_d ^ 
dt S^ 


+ 


■ Sf 


■ Sf 


+ 

V 5p\ 


Tp 


= ihi + (f) - H{z), 


SC ’ Sct^ ’ < 5 ^ 

such that the Euler-Poincare equations yield 

-Jz + JC - JC - {p\V^H{z)) = 0 

—ihp + —ihp] — [/?, ihC — H (z)] = 0 . 


M 

Sz 


= 0 . 


= -JC - (p|Vzi^(z)), 


thereby completeing the proof. ■ 

The observation that hybrid classical-quantum dynamics can be expressed by using the Heisenberg 
and unitary groups motivates us to investigate further the interplay between these two symmetry 
structures. The next Section shows that combining the two groups into a semidirect product yields 
the variational formulation of quantum expectation dynamics. 


4.2 The semidirect product (g) W(^) 

While the previous Section used the direct product x U{J^) group structure to obtain hy¬ 

brid classical-quantum dynamics, we shall now illustrate how constructing the semidirect product 
(S)7Y(J^) allows to shed new light on the dynamics of expectation values, thereby extending 
Ehrenfest theorem to more general situations. 

The semidirect product (S)iY(J^) can be constructed upon using the celebrated displace¬ 

ment operator from the theory of coherent quantum states. This is defined as follows 

Uhi’ix) = e-i V^(a; - h^), Vh = (h, ip) G , 

where the phase space vector h G is expressed as h = (hg,hp). This operator defines a group 
homomorphism that can be used to construct the following product rule in (^U(Jif): 

{hi,Ui){h2,U2) = {hh2,Ui{Uh,U2Ul)) , Vhi,/i2 G?^(m2-), V17i,C/ 2 GiY(jr), (51) 

where / 11/12 is the product rule in the Heisenberg group, already dehned in (48). Notice, upon denoting 
Z = {Q,P) (quantum canonical operators), the displacement operator Uh leads to the following Lie 
algebra homomorphism i : f)(M^"') —>■ u{Jif) 

c(C) = + c • J^), VC = (C, A € (i(m2-) , 

which occurs in the Lie bracket structure on (s)u(JP), given by 

^'^(Ci.a)(‘^2,C2) = (ad^iC2, [67(C2)] - [C27 (Ci)] + [6>6]) > 

where the operator ‘ad’ appearing in the first slot on the RHS is the infinitesimal adjoint action on 
f)(M^"'), as it was dehned in (50). No confusion should arise from this notation. 

In order to construct a dynamical theory by using the group structure above in the Lagrangian 
(45), we need to hnd an action of ?^(M^”)(DW(J^) on the space x S(JP). This task can be 
achieved by computing the coadjoint representation on the semidirect product. This computation can 
beneht from the following property. 
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Lemma 4 (Equivariance) With the notation above, the following relations hold 


UhZUl = Z-hi, i{Mh C) = Uh i{C) ul, 

where I is the identity operator on and Ad/iC = (C; i;^ + h-JI(^) is the adjoint representation on 


Proof. The first relation is easily proved by a direct verification. The first component reads as 
follows: 


{UhXUl)f;{x) = [Uhx) 








. hp cc 


2ft e " ft 'ijj{x + hq) 

. hp-(a;-hq) 


= e ^ e ^ 
= {x- hg)'ilj{x) 


{x - hg)e 


-I- 




Similarly, the second component reads 

{UhPul)if{x) = -ih{UhV) 


= -Uh 


e 

. hp • £C 


ft ll}{x + hq) 




eft'^e * 2ft e * ft (hp+ hg) + + h^)) 


= (-hp + P)V’(®) 


Combining both components proves the first relation in the lemma. The second relation follows by 
direct substitution 


t(Ad/, C) = i(c, 0 + h • JC) = ((j + h-K - Z -K) = -ifi i<p-iZ- hi) ■ JC) = 

=-ih[<f - {UhZuD ■K) = Uh{-ih{(j + iz-c))ul = Uhi{Oul, 


thereby completing the proof. ■ 

Eventually, by making use of the previous relations in the definition of coadjoint representation, one 
finds the following expression: 

Ad(;, f;)(z/,p) = (^- aSh + {fi- U^nU,ih-^SZ), a, ulu^ 

where we have used the notation n = {v,Oi) G ~ This coadjoint representation is 

computed explicilty in the Appendix A. 

Then, upon fixing the invariant set a = 1 and by introducing the variables z = —Jv and = 
ih~^p, we obtain the following action of on the space x 

<^ih,u){^,p^) = (z - h + {UZU^ - Z), ulu^p^UUh) , 

where we have used the expectation value notation (A) = {A\p^). 

4.3 Geometry of quantum expectation dynamics 

At this point, the semidirect product (^U{J^) has been characterized and it has been showed 

to possess an action on the classical-quantum phase space PJ^. Then, we consider the evolution 

of the classical-quantum variables {z, p^) under the action of {h~^,U~^), which then gives 

z{t) = zo + h(t) + {U{t)^ZU{t) - Z I p^g) , p^{t) = Uh{t)U{t)p^gU{tyUh{t)^ . (52) 
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The evolution above has the following crucial feature: 


z(t} - {Z\p^{t)) = zq- {Z\pp^) , 


as it is verified upon computing 

{Z I U{t)pp,U{t)^) = {Z I Uh{t)^p^{t)Uh{t)) = {Uh{t)ZUh{t)^ I p^{t)) = {Z- h{t)I I p^{t)). 

Therefore, in order to study expectation value dynamics, one can simply initiate the evolution under 
the initial condition zq = (Zlpp^), which is then replaced in ( 52 ). Moreover, the evolution above, 
produces the equations of motion 

^ = C - {[p,z\,C), p=[ih~^z-K + i^p] 

where C = h and ^ = llhUU^Ul- Analogous expressions hold for the variations {6z,6p). 

At this point, we consider the Euler-Poincare Lagrangian of the classical-quantum mean field 
model (49). Although that was written previously on the space (f)(M^"') x (M^”' x we 

now change perspective and we interpret the same expression (49) for p^) as a Lagrangian 

of the type 

(fi(M2")(Du(=^)) X X PJ^f) ^M. 

Notice that the Hamiltonian operator H{z) depends on the classical variable z, which has to be 
interpreted as the expectation value {Z\p^). This amounts to consider quantum systems for which the 
total energy can be written in terms of both the quantum state pp and its corresponding expectation 
values z = {Z\p^). (Notice that this is a very general case, as it shown by considering the kinetic 
energy expression ( P ^)^/2 = ( p )^/2 + {P - {p))^/2). 

Theorem 5 Consider the Lagrangian (49) and its associated variational principle for mixed quantum 
states ^ 

^ It + {pit)P^^it) +4>{t) --H'(z(i)))^dt = 0, 

with variations 

SC = 7, S(^ = d - c ■ Ij , 

SC = f/- ih~^ {[C,Z ■ J7] - [7?, Z • JC]) + [p, C ], 

= 7 - ([/9, Zfp) , 

5p = \ih~^Z • J7 + ry, p] , 

where 7, 0 and p are arbitrary and vanish at the endpoints. Then, this is equivalent to the following 
equations of motion 


z = JVz {p\H{z)) - ih-\z\ [H{z),p\ ) , (53) 

ihp = [H{z),p]+V^{p\H{z)) ■ [Z,p\ . (54) 
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Proof. 


This follows by a direct subsistution of the variations in the action principle. We have 


C{t) • Jz(t) + (p{t),ih^{t) + 4>{t) - H{z{t))) dt = 


j ^JIz • 5C + {p,S(j)) - {ihp,5^) - (JC + (/o|Vz^f(z))) • 6z + {ih^ + 0 - H{z),6p) ^ dt 
j (^Iz-j + (^p,9-C- J 7 ) - {ifip, V - ■ J 7 ] -[p,Z ■ JC]) + [??, ^]) + 

- (JC + {p\^xH{z))) ■ (7 - {[p,Z],p)) + {in^ + (t)- H{z), [ih~^ Z ■ Sj + 1], p])'^ dt 
I (( - Jz - {p\V,H{z)) + ( [ih-^p, H{z)] , JZ> , 7 ^ 

+ (i^P+ \p,Z ■ (p|Vzi?(z))j + [p,H{z)] ,r])'] dt 


Then, since 7 , 9 and p are arbitrary and vanish at the endpoints, the proof follows. ■ 

In order to understand how the above result is related to the usual Ehrenfest equations for quantum 
expectation dynamics, we immediately observe how these equations (53)-(54) are recovered (along 
with the evolution of p) in the case when Vz^I(z) = 0. As it was pointed out previously, the new 
feature of equations (53)-(54) lies in the fact that the expectation values have been considered as 
independent variables already occurring in the expression of the conserved total energy {H{z)). This 
confers the system (53)-(54) a hybrid classical-quantum structure. Indeed, one observes that new 
coupled classical-quantum terms appear in Ehrenfest dynamics: these are the first term on the RHS 
of (53) and the second term on the RHS of (54). 

Notice, the first term on the RHS of (53) does not involve the quantum scales given by h. Eor 
example, a purely classical system is given by a quantum phase-type Hamiltonian operator of the form 
H[z) = /i(z)l, where h{z) is the classical expression of the Hamiltonian. In this case, while equation 
(53) recovers classical Hamilton’s equations, the quantum evolution (54) specializes to coherent state 
dynamics of the type 

ihp = Vxh ■ [Z, p] . 

This establishes how quantum states evolve under the action of purely classical degrees of freedom, 
thereby enlightening the interplay between classical and quantum dynamics. The same equation can 
also be obtained by linearizing the quantum Hamiltonian operator H(Z) around the expectation values 
(i.e. in the limit Z —)■ zl), as prescribed by Littlejohn’s nearby orbit approximation for semiclassical 
mechanics [29]. 


5 Conclusions 

This paper investigated the geometric symmetry properties of quantum and classical-quantum vari¬ 
ational principles. Upon departing from the Dirac-Erenkel Lagrangian, different quantum mechanics 
pictures were recovered from the same variational principle, upon making extensive use of Euler- 
Poincare theory. This reduction by symmetry naturally leads to consider the Hopf bundle as the 
natural setting for the Schrodinger picture of pure state dynamics, as already proposed by Kibble 
[27]. In addition, new variational principles were presented for mixed state dynamics in both the 
density matrix and the Wigner-Moyal formulation. Later, new quantum variational principles were 
also presented for the Heisenberg and Dirac’s interaction pictures of quantum dynamics, where the 
isotropies of the initial state was shown to possess the same role as phases in the Schrodinger picture. 
In particular, Dirac’s interaction picture involves the geometric semidirect product structure of two 
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different unitary groups associated to the different quantum propagators arising from the splitting of 
the Hamiltonian operator. 

In the last part of the paper, the Dirac-Frenkel Lagrangian was augmented to account for classical 
degrees of freedom, incorporating the dynamics of classical motion in hybrid classical-quantum dy¬ 
namics. As a first step, the mean field model of classical-quantum dynamics was described by using the 
direct product of the Heisenberg group (governing classical evolution) and the unitary group 

(governing quantum evolution). Later, Ehrenfest’s expectation value dynamics was shown to 
arise from a novel set of equations, whose dynamics evolves both expectation values and the quantum 
density matrix under the action of the semidirect product ?^(M^”’)(S) whose group and Lie 

algebra structures were presented. In this context, purely classical dynamics was shown to arise by 
using Littlejohn’s nearby orbit approximation [29]. 
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A Adjoint and coadjoint representations of 


First, by using the product rule (51), one computes the explicit formula for the conjugation action of 
(J^) on itself 

\h,,U,){h2,U2) = {hi,Ui){h2,U2)ihi,Ui)-^ = 

= ((hi,^i),t/i)((h2,992),t/2)((-hi,-99i),[/t^l7|l7fc,) = 

= ((hi,^i),t/i)((h2 - ih2 • JIhi),D2 (Uh,{UlulUh,)Ul) ) = 

= ((h2,(/P2 -h2 • Jhi), UiUhMUh.ululUh.uljul ). 

Then, taking an arbitrary curve 


{h 2 {t),ip 2 {t),U 2 {t)) such that (h2(0), (^2(0), D2(0)) = (0,0,1) , 

and upon denoting (h2(0), </32(0), 1/2(0 )) = (Cj </’)?) ^ fl(M^"') (Du(J^) and Uh 2 { 0 ) = i{C), one defines 
the adjoint action of on its Lie algebra as follows 




_d 

dt 

_d 

dt 


t=o 


I(/ii,i7i)(^2(l), 1 / 2 ( 1 )) - 


t=0 


((h2(l),(^2(l) -h2(l) -Jhi), UiUh,U2{t){UH2{t)ululUh,ulii))ul 

h2{0),M0) - MO) ■ Sh,), 

+ 1/1 Uh, U2 (0) {Uh2 ml ul Uh, ul (0)) ul + 

- uiumo){UhMuluMMh2mh,mimui) = 


((c,<^ - c • jhi), umc+muiui - uhaouI) = 

(^AdhC, uium+muiui-tiAdhO) • 
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At this point, using the notation (z/, /i) = x ~ x u(Jf')*, one 

computes the coadjoint representation on via the pairing 

( (C>0) = Ad(;,_[/)(C,0) = 

= ((1.,/i), (Ad,c, UiUh,{i + i{Q))ulul - UhAOul)) = 

= ((Ad;; z. + (ul{u^^lU - fi)Uh ), ulu^fiu^ut), (c,o) , 

where Ad;; z^ = {u — ajh, a) is the coadjoint representation on and l* : is 

the dual of the Lie algebra homomorphism generated by the displacement operator Uh, which is also 
computed via the pairing giving the following expression 

= (^(/r,-z^-ijZ),Tr(z^-V)) , G 

By direct substitution, the coadjoint action on (s)ZL(^) reads 

Adli^ u^n) = aJh + {n- U^nU,ih-^SZ), a, ulu^fiUUh^ . 
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